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Because of its relatively simple structure the Nambu-Jona-Lasinio (NJL) model p] is 
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We study the conditions for the existence of stable quark matter in the Nambu-Jona- 
Lasinio mean field at zero temperature and discuss its interpretation. 
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(N . 1. INTRODUCTION 

<N ■ 
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one of the most popular models for studying the spontaneous breaking of chiral symmetry 
and its restoration at finite temperatures or densities [[J - @. However, most of its 
predictions suffer from the fact that the model is a pure quark model without confinement, 
which in general leads to rather unrealistic scenarios. In mean field approximation for 
q^. T = we can distinguish three different cases 0: 

Case I: The phase transition is of second order. Then for all densities there is a uniform 
phase of quarks. The quarks are massive at low densities and (almost) massless at high 
densities. The pressure is always positive and if we do not apply an external force to the 
^ ■ system it will expand to arbitrarily large volumes, i.e. arbitrarily small densities. 

Case II: The phase transition is of first order. In this case there will be a regime in 
which two phases coexist. In general, this mixed phase consists of a low-density phase 
(p = pi) with massive and a high density phase (p = ph) with massless quarks. However, 
for p < pi we have again a uniform phase of massive quarks. As in case I, low densities 
are energetically favored and quark matter is unstable against expansion. 

Case III: The phase transition is of first order, but with pi = 0. Then the low-density 
component of the mixed phase is the vacuum, i.e. we find droplets of massless quarks 
surrounded by the non-trivial vacuum. At least in a schematic way this is very reminiscent 
of the MIT bag model. The droplets are the energetically most favored configuration and 
therefore stable against expansion or collapse. In contrast, any uniform distribution of 
quarks at lower densities is unstable against break-up into the droplet-vacuum phase. 
This scenario was also studied in ref. || within a similar model. 

For NJL-like models it depends on the model parameters which of the three cases is 
realized. However, the non-existence of a uniform dilute gas of quarks in nature excludes 
the cases I and II, which both predict such a phase, from being realistic. In case III there 
is no stable uniform quark gas at low densities but instead we have the droplet phase. At 
least if we adopt the interpretation of quark droplets as schematic bag model baryons this 
seems to be the most realistic scenario of the model. In this article we want to discuss 
whether it is also compatible with "realistic" model parameters. To large extent this will 
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be done following ref. 0. 

2. FORMALISM 

We consider a Lagrangian for quarks with rif = 2 flavors and n c = 3 colors interacting 
by NJL-like four-fermion vertices: 

£ = $(tf-mo)i> + G S [(W) 2 + $ij 5 rijj) 2 ] - GV(WV) 2 • (1) 

The first two interaction terms (scalar-isoscalar and pseudoscalar-isovector) were taken 
from the original N JL-Lagrangian [|l| . Since vector interactions are known to be important 
at finite densities, e.g. like in the Walecka model |5J, we have also included a vector- 
isoscalar term. Linearising t/jip and if)'ij) about their thermal expectation values we can 
calculate the mean field thermodynamic potential (per volume) at temperature T and 
chemical potential fi Q. We restrict ourselves to the Hartree approximation. The result 
for T = and \i > reads: 

u MF (p; m^ R ) = -12 / jZ^fa+fa-EjOfa-Ep)] + {m ~™f ~ {jL ^r^ , (2) 

with E p = \Jrri 1 + p 2 . The integral is strongly divergent and has to be regularized. For 
simplicity we use a 3-dimensional sharp cutoff. The auxiliary variables m ("constituent 
quark mass") and hr are defined as m = m — 2Gs{ipip) and //r = h — 2Gy{i()ty). Since 
on the other hand the condensates (ipip) and {ip^fp) have to be calculated from the ther- 
modynamic potential by taking the appropriate derivatives we encounter a selfconsistency 
problem. It can be shown that it is solved by the stationary points of ujmf with respect to 
m and fj,R. This can be used to eliminate hr for given values of \i and m, fiR = iar(h, m). 
Furthermore, we subtract a constant in order to shift the value of the minimum of the 
vacuum thermodynamic potential to zero, i.e. we define: 

u)(fi,m) := wmf(k m,/i R (/x,m)) - uj M f(0] m vac ,0) . (3) 

Here m vac is the constituent quark mass which minimizes the thermodynamic potential 
in vacuum. For a given value of /x the extrema of Co correspond to the selfconsistent 
solutions and the stable solution is given by the absolute minimum. Other thermodynamic 
quantities, like the pressure p, baryon number density Pb = ^- (^V) an d energy density 
e, can be calculated from the thermodynamic potential in the standard way. 

3. STABILITY OF QUARK MATTER AT T=0 

With the appropriate choice of parameters chiral symmetry is spontaneously broken in 
the NJL-vacuum, i.e. the quarks acquire a constituent mass m vac much larger than the 
current mass mo- Unless stated otherwise we will work in the chiral limit, mo = 0. Then 
the thermodynamic potential u is symmetric about m = 0. In vacuum m = corresponds 
to a local maximum and Co is minimal at m — ±m mc . For large chemical potentials chiral 
symmetry becomes restored. Thus, with increasing /i, the maximum at m = eventually 
has to turn into an absolute minimum. If this happens at a critical chemical potential less 
than the vacuum mass, p crit < m vac , we have a first-order phase transition of the type 
"case III". 
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An example for this case is shown in fig. [I] where we find a critical chemical potential 
f^crit = 376 MeV, well below the vacuum mass m vac = 400 MeV. Since for T = and 
Gy = the chemical potential is equal to the Fermi energy of the system, for a given 
mass m, a chemical potential less than m corresponds to zero density and vice versa. 
Hence the degenerate minima of u at fi — p cr u correspond to two phases of equal pressure 
in chemical equilibrium, one being the nontrivial vacuum (m = m vac and p# = 0) and 
the second consisting of massless quarks at finite density. In our example we find p# = 
0.47/m~ 3 = 2.75po for the dense phase. This means, droplets of massless quarks and a 
baryon number density of 2.75po are stable in vacuum, as far as our mean field approach 
is justified. 




Figure 1. Thermodynamic potential wasa 
function of the constituent quark mass m 
for m = 0, A = 587.9 MeV, G S A 2 = 2.48 
and Gy = 0. The different curves cor- 
respond to three different chemical poten- 
tials: /i = (dashed line), p = 376 MeV 
(solid line) and p = 450 MeV (dotted line). 



Figure 2. Energy per baryon number 
function of the baryon number density ps 
for the same parameters as in fig. |l|. The 
solid line corresponds to massive, the dot- 
ted line to massless quarks. 



Alternatively this can be seen if we calculate the energy per baryon number, e/pb as a 
function of density ps- The result is shown in fig. |2| for the same set of parameters as in 
fig. III. At not too high densities there are solutions with massive and with massless quarks, 
and the massive ones are always energetically favored. In the limit pb — > the energy per 
baryon number is just 3 times the constituent quark mass. However, in agreement with 
our considerations above we find the absolute minimum of e/ps at 2.75po which is in the 
regime where only massless solutions exist. At this point the pressure vanishes and the 
energy per baryon number is 3 times p cr it- This shows again that p crit < m vac is necessary 
for this minimum to be stable. 
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The results can become qualitatively different if we change the model parameters. For 
instance, with the parameters of fig. [I] but with a vector coupling Gy = 0.5Gs we find a 
first-order phase transition at fi cr i t = 410.3 MeV > m vac . This corresponds to what we 
called "case II" and there is no dense matter solution which can coexist with the vacuum. 
The energy per baryon number as a function of density looks similar to fig. 0, but with 
the minimum of the massless solutions being only metastable: The corresponding energy, 
3/i cr j t , is now larger than the energy per baryon number for p B — > 0, which is 3m vac . 
Finally, if we further increase the vector coupling, e.g. Gy = Gs, we find a second-order 
phase transition and the energy per baryon number is a strictly rising function of density. 

A more systematic overview is given in fig. ||| where three lines of constant binding 
energy per quark, Eu n d = mvac — parity are plotted. For each value of m vac the cutoff A 
was fixed by fitting the pion decay constant to its empirical value f n = 92.4 MeV. We 
are left with two model parameters, which can be chosen to be m vac and the ratio of 
vector and scalar coupling constant, Gy/G$- We see that the vector interaction reduces 
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Figure 3. Lines of constant binding energy 
per quark for fixed f n = 92.4 MeV and 
varying constituent quark masses m vac and 
coupling constants Gy/Gs- MeV (solid 
line), 50 MeV (dashed) and 100 MeV (dot- 
ted). 



Figure 4. Energy per baryon number as 
a function of baryon number density for 
non-strange quark matter in the 3-flavor 
model. The parameters are: = mjj = 
5.5 MeV, m s = 140.7 MeV, A = 602.3 MeV 
and G S A 2 = 1.84, KA 5 = 12 A (solid line), 
G S A. 2 = 2.31, K = (dashed line). 



the binding, which is not surprising since vector mean fields are known to be repulsive. 
We also find that the binding becomes stronger with increasing constituent quark masses. 
For m vac < 343 MeV there is no bound quark matter, even without vector interaction. 
This agrees quite well with an approximate relation derived in ref. [0] which states that 
stable quark matter is only possible in the NJL-model if m vac > 4f n . 



5 



The region shown in fig. ^| corresponds more or less to the regime of "realistic parame- 
ters" . One could try to further constrain m vac by fitting the quark condensate. However, 
{ijjip) is not known very precisely and its dependence on m vac is only weak, once f n is 
fixed. The vector coupling constant could be determined by fitting vector meson masses. 
In the literature this usually leads to Gy/Gs of the order of 0.5 to 1 @|| although much 
higher values are also found |TD|. On the other hand a realistic vector coupling constant 



in dense matter could be rather different from the vacuum one. 

Thus the scenario of stable quark droplets in vacuum ("case III") is possible but not 
a necessary or very probable consequence of restricting the model parameters to realistic 
values. This is in some contradiction to ref. || where stable droplets are found to be 
most likely. The model of ref. M is almost identical to our model with Gy = 0, with the 
only difference that instead of using a sharp cutoff the vertices are multiplied with form 
factors of the form ( p 2 A [ _ A 2 )" • ^ s an exam pl e we looked at v = 1 and A = 800 MeV. It 
turned out that we are able to produce very similar results with our model if we use a 
sharp cutoff of about 350 MeV. However, for this cutoff parameter we find a pion decay 
constant of less than 60 MeV. This suggests that the parameters chosen in ref. || might 
also correspond to a too small value of /„-. 

4. QUARK DROPLETS AS SCHEMATIC MIT-BAGS 

In case III dense matter of massless quarks can coexist with the non-trivial vacuum. 
As pointed out in the introduction it is very attractive to identify these quark droplets 
with baryons in a bag-model picture. In fact, the energy of a spherical MIT-bag [ITT] with 
radius R, 

47T 3X - Z B , n . /4tT\ 1/3 1/3 

Emit = -j R B + ~ ^ = ~ + ( 3x " *o) ( y j Pb , (4) 

has a similar structure as the energy per baryon number for massless quarks in the NJL 
mean field 



2 \ 1/3 

Ps /3 + G v n 2 c p B , (5) 



e B 3n c I 3tt 
Pb Pb 4 ' 

if we switch off the vector interaction. Here the bag constant is given by B = £(0,0). 
Of course our thermodynamic approach is valid only for a large number of particles in 
a large volume whereas eq. (|]) was derived for three quarks in a volume which can be 

1 /3 

small. This leads to different coefficients in front of p B . However, taking x = 2.04 and 



Zq = 1.84 [iTJ, the deviation turns out to be surprisingly small (~ 20%). One should 
also keep in mind that eq. ([5]) is only valid for the massless quark solutions. As we have 
seen, if the density is low enough the quarks can lower their energy by acquiring a finite 
constituent mass (cf. fig. |^). In particular, only a finite amount of energy is needed to 
lower the density to zero, corresponding to an infinite bag radius: the NJL mean field 
does not confine. As an important consequence the quark droplets become immediately 
unstable against evaporation of massive quarks at any finite temperature and one has to 
be very careful with the interpretation of the model for T > 0. 
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5. THREE-FLAVOR MODEL 



In this section we extend the model to three quark flavors 12]: 



k=0 

- K[det f $(l + l5 )i>)+det f $(l - 7b )V)] • (6) 

Here ip = (u, d, s) T is a 3- dimensional vector and m = diag(m^, Wq, tuq) a 3 x 3 matrix in 
flavor space. In addition to the 4-point interaction the Lagrangian contains a t'Hooft-type 
6-point interaction which is a determinant in flavor space and breaks the Ua(1) symmetry. 
For simplicity we neglect vector interactions in this section. Since each flavor is separately 
conserved we have in principle three different chemical potentials and the thermodynamic 
potential contains the term fj,N = p u u^u + fidd'd + psS^s. [] 

In the following we consider non-strange isospin-symmetric matter, mj} = ttIq, p u = pd 
and p s — 0. We adopt the model parameters of ref. [13] which have been fitted to the 



pseudoscalar mass spectrum and f v . For these parameters we find stable quark matter 
with p = 2.36po an d a binding energy per quark of 6.5 MeV. The energy per baryon 
number is plotted in fig. |4] (solid line). Because of the determinant interaction this result 
is influenced by the strange quark even though the net strangeness (s's) is zero. To show 
this we switch off the determinant interaction and increase G$ such that the vacuum 
properties of the non-strange sector, (mm), f n and m n , remain unchanged. The result is 
also shown in fig. |] (dashed line). Here the binding energy is 2.9 MeV at p = 2.43po- 
Thus the flavor mixing gives some extra binding but the effect is not very large. A more 
systematic discussion of the SU(3) model will be published elsewhere. 
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1 Here we use the notation of ref. ||. Notice that /i s is different from the strangeness chemical potential 
Estrange which is often introduced as the deviation from a common baryon chemical potential fiB- I^N = 
Hb{u u + d'd + s's) + UstrangeS^ s. In particular in our notation /j, s = corresponds to (s^s) = 0. 



